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The Tolman VII solution is an exact static spherically symmetric perfect fluid solution of Ein- 
stein's equations that exhibits a surprisingly good approximation to a neutron star. We show that 
this solution exhibits trapped null orbits in a causal region even for a tenuity (total radius to mass 
ratio) > 3. In this region the dynamical part of the potential for axial w - modes dominates over 
the centrifugal part. 
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Exact solutions of Einstein's equations provide both a route to the physical understanding of relativistic phenomena 
and a check on numerical procedures. Unfortunately, even in the case of isolated static spherically symmetric perfect 
fluids, very few solutions are known which pass even elementary tests of physical relevance [Q. In a recent study 
H-j | H four types of ostensibly reasonable solutions were examined for the internal trapping of null geodesies along 
■ with the associated resonance scattering of axial gravitational waves Q , || , || . The solutions examined iri|| have 
a large abrupt density discontinuity ("phase transition") at the boundary 0. As a consequence, though physically 
• reasonable, they are not good models of (say) neutron star structure There are only two known exact static 
spherically symmetric perfect fluid solutions which match onto vacuum without any density discontinuity at the 
boundary: the Buchdahl solution JoJ , and a special case of the Tolman VII solution j[o|. These solutions are both 
quite reasonable models of neutron star structure j|] . The Buchdahl solution does not exhibit internal trapping of null 
geodesies || . In this paper we show that the Tolman VII solution exhibits trapped null orbits even for a tenuity (total 
radius to mass ratio) > 3. In this region the dynamical part of the potential for axial w - modes actually dominates 
y—i ' over the centrifugal part. 

o : 
o ' 

qh. II. THE TOLMAN VII SOLUTION 

6J[). A. Outline of the solution 

> 

X 

ds 2 = ds 2 T +r 2 dn 2 , (1) 

where dfl 2 is the metric of a unit sphere and df, is a static Lorentzian 2-surface with coordinates {t,r). The source of 
(0), by way of Einstein's equations, is considered to be a perfect mathematical fluid generated by radial streamlines 
of constant r. The solution terminates and matches onto a Schwarzschild vacuum upon the vanishing of the isotropic 
pressure. 

The Tolman VII solution begins with the ansatz 

8n P^ - 1 (3 2 T - (-) 2 ) (2) 

where p(r) is the energy density and /3, T and K are constants. The physical meaning of these constants is clear at 
the outset: K 2 — 8 J* , where p is the central density, and p = jgij for = (3. The isotropic pressure (p) vanishes 
at = (3 defining the boundary (£). (The Schwarzschild interior solution is not the (3 — > oo limit.) Physically 
reasonable values of T for neutron stars are T = 1 + e where e <~ 1CP 8 . In this paper, without any change to our 
conclusions, we take T = 1. It is convenient to use the dimensionlcss variable x = 
With (J2j) the effective gravitational mass (to) is given by 

m(P, x ) = K^(5-^x 2 ) (3) 



We begin with the static spherically symmetric metric in "curvature" coordinates |TT 
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which gives a tenuity ^j7^y|s = -ps = ol. In contrast to @ and (|J), the isotropic pressure is not a simple function. 
(Indeed, this caused Tolman to abandon the solution [fo|.) The pressure is given in the Appendix. We record the 
metric: 

d Sr = -D „„(I n (y ^ )) df + w - y (4) 

where 

a{[3,x) = l-x 2 + —, (5) 

^)^/f(* 2 -^), (6) 

and D defines the choice of scale for t (which, of course, is of no physical consequence). The number C(/3) (which is 
governed by T) is given in the Appendix. Without regard to causality, the allowed range in (3 is < (3 <~ 1.3896 
where the upper bound is set by the requirement of a finite central pressure. 



B. Notation 



We simply review here some key equations to set the notation (see e.g. [|| for details). Non-radial null geodesies 
satisfy 



rV* = (l-^)(«-l), 



with 



B(r) = ~ (8) 
V9tt 

where • = d/dX for affine A, and 6 is a constant > 0, the "impact parameter". The evolution of null geodesies is 
restricted by the condition b < B(r). The odd parity (axial) w-modes are non-radial perturbations of the spacetime 
which do not couple to the fluid at all. In terms of the frequency w and mode number I > 2 the governing equation 
is given by jj) 

(J^+^Z^V^Z, (9) 

where r* is the "tortoise" coordinate. The potential is conveniently expressed in terms of r and is given by 

V M = 1^2 W + !) + 4 ^ 2 (PM - p(r)) - ^). (10) 

The "centrifugal" part of the potential is given by ^fe^ , the remainder being the "dynamical" part. A necessary 
condition for the occurrence of resonance scattering of axial gravitational waves by an isolated distribution of fluid is 
a local minimum in V(r) within the boundary of the fluid. 



C. Properties 

We start by noting that for an appropriate choice of the constant D in metric (||) the potential impact parameter B{r) 
and associated potential V(r) match continuously onto their exterior Schwarzschild vacuum values at the boundary S. 
The appropriate value of D is shown in Figure 1 where, for convenience, we also record the tenuity m ^ r ^ |g = a. The 
choice for D is of no physical consequence. We have included Figure 1 simply to emphasize the fact that the continuous 
matching of B(r) and V(r) is possible. Figure 2 shows the value of ™ for which B(r) reaches a local maximum 
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compared to the value of below which the classical adiabatic sound speed (v^ = becomes superluminal. We 
conclude that (3 <~ 1.29663 and that the tenuity >~ 2.97402 for the trapping of null geodesies in a causal (v s < 1) 
region. Figure 3 records some gross features of the potential impact parameter B(r). The minimum (3 for which 
there is a local maximum in B{r) is (3 ~ 1.24955 corresponding to a maximum tenuity of ~ 3.20229. The Tolman 
VII solution therefore exhibits trapped null orbits for a tenuity > 3. This is the first exact causal solution which 
has been shown to have this property. As regards the full potential V(r), we find that the minimum (3 for which 
there is a local minimum in V(r) is /3 ~ 1.28582 corresponding to a maximum tenuity of ~ 3.02419. This means 
that the dynamical part of the potential for axial w - modes in fact dominates over the centrifugal part in the range 
of tenuity ~ 3.02419 < a <^ 3.20229. This is examined in Figure 4. Some gross features of the potential V(r) are 
shown in Figure 5. For completeness, a comparison of the Tolman VII solution with the well known (and widely used) 
Schwarzschild interior solution is given in Figure 6. 

III. DISCUSSION 

The Tolman VII solution is the first exact causal solution which has been shown to exhibit trapped null orbits for 
a tenuity > 3. Moreover, this solution shows that the dynamical part of the potential for axial w - modes can in fact 
dominate over the centrifugal part over a range in tenuity. These features would not be expected from any studies of 
the Schwarzschild interior solution. Whereas the potential V(r) satisfies the necessary condition for the existence of 
trapped w-modes even for a tenuity > 3, it seems unlikely that trapped w-modes actually exist in this solution for 
tenuities in this range. Nonetheless, this solution makes it clear that trapped null orbits and w-modes exist in causal 
solutions with realistic equations of state and everywhere smooth densities. There is reason to believe then that the 
purely relativistic phenomena of trapped null orbits and w-modes are part of real astrophysical systems and not just 
curiosities of toy models. 

ACKNOWLEDGMENTS 

This work was supported by a grant (to KL) from the Natural Sciences and Engineering Research Council of 
Canada. Portions of this work were made possible by use of GRTensorll |fl2|| . We would like to thank Don Page for 
several very helpful comments on the Tolman VII solution. 



Electronic Address: lake@astro.queensu.ca 



[1] M. Delgaty and K. Lake, Comput. Phys. Commun. 115, 395 (1998) ( |gr-qc/9809013[ ). This work (with corrections and 

updates) is being incorporated into an on-line database ( http: //grtensor .phy . queensu. ca/ so lutions/ ). 
[2] M. Ishak, L. Chamandy, N. Neary and K. Lake, Phys. Rev. D 64 024005 (2001) ( gr-qc/00070"73| ). 

[3] The trapping of null geodesies inside the Schwarzschild interior solution was studied long ago. (See, for example, F. de Felice, 
Nuovo Cimento 63B, 649 (1969) and J. Jaffe, Mon. Not. R. Astron. Soc. 149, 395 (1970)). The trapping of null geodesies 
requires a detailed consideration of scattering. 

[4] S. Chandrasekhar and V. Ferrari, Proc. R. Soc. Lond. A 434, 449 (1991). See also N. Andersson, Y. Kojima 
and K. D. Kokkotas, Ap. J. 462, 855 (1996) and N. Andersson and K. D. Kokkotas, Mon. Not. R. As- 
tron. Soc. 297, 493 (1998). An extensive review relevant to the background of this paper is given by 
K. D. Kokkotas and B. G. Schmidt "Quasi-Normal Modes of Stars and Black Holes", Living Reviews 1999-2 
(http : //www. livingreviews . org/Articles/Volume2/1999-2kokkotas/ ). 

[5] The idea of trapped waves in compact objects is not new. See B. Iyer and C. Vishveshwara in A random walk in relativity 
and cosmology Ed: N. Dadhich, J. Krishna Rao, J. Narlikar and C. Vishveshwara, (Wiley Eastern Ltd, New Delhi, 1985) 
and references therein. 

[6] It is now well known that sufficiently compact polytropic and uniform density stars can support " w-modes". The polytropic 



case is essentially numerical in nature. See m and K. Rosquist, Phys. Rev. D 59 044022 (1999) (gr-qc/9809033). The case 



examined by Rosquist is acausal in the regi on of interest. A lthough there is much work on the relatively simple uniform 



density case (See, for example, Z. Andrade (jr-qc/0103062) and references therein), the model is but a "toy" due to the 



everywhere infinite sound speed. Similarly, the Zel'dovich "stiff matter" models (with everywhere lightlike sound speed) can 



3 



be used only as a guide to demonstrate what might be possible in reality. See M. Karlovini, K. Rosquist and L. Samuelsson, 
Class. Quant. Grav. 18, 817 (2001) and references therein. (The authors construct "toy" models with multiple necks.) 

[7] Whereas the density a t the surface i s not required to be (see, for example, P. Musgrave and K. Lake, Class. Quant. 
Grav. 13, 1885 (1996) ( jr-qc/9510052)) models which are otherwise physically acceptable may give physically unreasonable 
values of the surface density. Uniform density models are, in this sense, the worst possible models. The surface to central 
densities for the models studied in || are: 0.17, 0.38 and 0.46 for the Finch-Skea, Heintzman and Durgapal solutions 
respectively. These values are about seven ord ers of magnitude t oo large for neutron stars. 

[8] J. Lattimer and M. Prakash, Ap.J. (in press) ( |astro-ph/0002232[ ). 

[9] H. A. Buchdahl, Ap.J., 147, 310 (1967). 
[10] R. C. Tolman, Phys. Rev., 55, 364 (1939). Mehra's solution (A. L. Mehra, J. Aust. Math. Soc. 6, 153 (1966)) is a special 

case of the Tolman VII solution. 
[11] We use geometrical units throughout. For example, the tenuity r-s/M is actually c 2 rs/GM. 

[12] Portions of this work were made possible by use of GRTensorll. This is a package which runs within Maple. It is en- 
tirely distinct from packages distributed with Maple and must be obtained in dependently. The GRTensorll softw are and 
documentation is distributed freely on the World- Wide- Web from the address tittp: //grtensor .phy .queensu. ca 



APPENDIX A: ISOTROPIC PRESSURE 



The pressure can be given in the form 



1 2xVa(^) /3 y/a(fcx)+b(P,x) , 
8np{p,x) = ( ^-cot(lnd — )+a{fi,x)-l), (Al) 

where a((3,x) and b((3,x) are given by (^]) and (^J) respectively and C(/3) is a constant (determined by the condition 
T = 1) given by 

C((3) = -^=exp(2arctan^(f3))(S + f ), (A2) 

where 



and 



Note that p{x =(3)=0. 



6 = y/15 - 6/3 2 (A3) 



= 23/3^ + 24/3*- 60(0 + g) 
[P) ~ /3(-23/3 2 + 60 + 4f35) ' [ 1 
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FIG. 1. The abscissa is (3. The upper curve gives the value of D 2 for which the potential impact parameter B(r) and 
associated potential for axial w - modes V(r) match continuously onto their exterior Schwarzschild vacuum values at the 
boundary E (r = f3K). The scale for D 2 is on the left. The lower curve gives the tenuity ^tj|e = a for which the scale is on 
the right. 

FIG. 2. The ordinate is and the abscissa (3. The curve uppermost on the left gives the value of for which the potential 
impact parameter B(r) reaches a local maximum. The other curve gives the value of below which the sound speed becomes 
superluminal. The curves cross at (3 ~ 1.29663 at a tenuity ~ 2.97402 

FIG. 3. Gross features of the potential impact parameter B(r). The vertical scale is irrelevant. The abscissa is -Jj^. The 
minimum f3 for which there is a local maximum in B(r) is j3 ~ 1.24955 corresponding to a maximal tenuity of ~ 3.20229. (a) The 
values of (3 and tenuity top to bottom are (1.29663, 2.97402), (1.29, 3.005), (1.28, 3.05176), (1.27, 3.1), (1.26, 3.14941), (1.25, 3.2). 
(b) An enlarged view of the lowest curve in (a). 

FIG. 4. Comparison of the potential impact parameter B(r) and the potential V(r). The abscissa is The vertical scale 
is irrelevant, (a) Dominance of the dynamical part of the potential. Here (3 = 1.28. The potential V(r) (upper curve on the 
left) is monotone decreasing yet the potential impact parameter B(r) reaches a local maximum in a causal region, (b) Here 
f3 — 1.29663 and B(r) (upper curve on the left) has a local maximum where the sound speed becomes luminal. 

FIG. 5. Gross features of the potential V(r). The abscissa is The vertical scale is irrelevant. The minimum j3 

for which there is a local minimum in V(r) is j3 ~ 1.28582 corresponding to a maximal tenuity of ~ 3.02419. (a) The 
values of (3 and tenuity top to bottom are (1.289, 3.00929), (1.29, 3.00463), (1.291, 2.99997), (1.292, 2.99533), (1.293, 2.99070), 
(1.294, 2.98608), (1.295, 2.98147), (1.296, 2.97687), (1.2966, 2.97415). (b) An enlarged view of the uppermost curve in (a). 

FIG. 6. The Tolman VII solution compared to the Schwarzschild interior solution. In each case the abscissa is a dimensionless 
measure of distance to the boundary E (set at 1). The vertical scale is irrelevant. The sound velocity is subluminal throughout 
the range shown for the Tolman VII solution (and infinite for the Schwarzschild interior solution), (a) The potential impact 
parameter B(r) for the Tolman VII solution (upper curve) compared to the Schwarzschild interior solution. The tenuity for 
the Tolman VII solution is 3.052 compared to 2.632 for the Schwarzschild solution, (b) The potential V(r) for the Tolman 
VII solution (lower curve) compared to the Schwarzschild interior solution. The tenuity for the Tolman VII solution is 2.974 
compared to 2.632 for the Schwarzschild solution. 
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